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Abstract. The paper deals with "quantization" and "second quantization" of 
phenomenological thermodynamics with respect to the Boltzmann's constant. It 
is suggested to perceive the quasithcrmodynamic parameter (corresponding to the 
Boltzmann's constant) as a mathematical analogue of the semiclassical param- 
eter (corresponding to the Planck's constant), and to introduce a new concept 
of a "thermocorpuscle" (a thermodynamic analogue of a particle where the co- 
ordinates are replaced by the nonequilibrium thermodynamic forces and the mo- 
menta are replaced by the corresponding flows). The semiclassical quantization of 
phenomenological thermodynamic Lagrangian manifolds yields a new system of 
equations for the quasithermodynamic fluctuations along a curve of evolution of 
a nonequilibrium physical system. This leads to a quasithermodynamic analogue 
of Bell's inequalities and their violation is a new effect that can be tested exper- 
imentally. The generating function of the quasithermodynamic fluctuations (the 
nonequilibrium analogue of a partition function) is interpreted as an expectation 
value of a second quantized operator expressed via the density of "thermocorpus- 
cles" . An analogue of the BBGKY chain of equations defines a deformation of 
the fluctuations by an interaction between the thermocorpuscles. In place of an 
interaction parameter in mechanics (the "external" Planck's constant), one intro- 
duces the "external" Boltzmann's constant for an asymptotic expansion of the 
thermodynamic collision integral. 



1. Introduction 

The original motivation for this paper stems from an idea to study the quasither- 
modynamic fluctuations in nonequilibrium statistical physics using the methods of 
semiclassical approximation of quantum theory. If we look at a semiclassical wave 
function i/jh(x), of a quantum system with n degrees of freedom 
and take a classical limit H — > for the corresponding Wigner's quasiprobability 
function p[ipf^(x,p), p = (pi,P2, ■ ■ ■ ,Pn), then under some mild conditions and as- 
sumptions [Bj, we obtain a smooth manifold 

A := supp \im p[ipti](x,p), 
h— >o 

where the limit on the right-hand side is understood in the weak sense. This manifold 
is a Lagrangian submanifold A C K^™ in the classical phase space with respect 
to the canonical symplectic structure co = Yli=i A dxi. 

The Lagrangian manifolds arise in a perfectly natural way in the equilibrium 
thermodynamics as well. Consider, for example, a physical system consisting of 
v moles of a chemical substance in a closed volume V. Then what happens in the 
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molecular kinetic theory is that one splits v into a product of a "very small" quantity 
k B (the Boltzmann's constant) and a "very big" quantity N (the number of particles 
in the system): 

v = ^k B N, (1) 

where R is a constant fixing the units of measurement (the universal gas constant). 
The basic idea is to consider a family of underlying systems with N\ = XN particles 
and volume V\ = XV, where A is the rescaling parameter, and to derive the em- 
pirical laws from the thermodynamic limit X — > oo. As a model example, take the 
underlying system to be a system of N\ particles on a n-dimensional torus of radius 
L\ = X l l n L, V = L n , L is fixed. Assume that it is described by an energy spectrum 
< E± < - • • < Em < . . . , where each level Em\ m 6 Z^ , has a finite degree 
of degeneracy gin* ■ One can construct a partition function (assume that the power 
series converges): 



oo _ 



m=0 



where j3 = T 1 is the inverse absolute temperature of the system. The free energy 
JF^{T) of the system at temperature T is given by 

4 X J(T) :=-k B T\nZ^(T- 1 ). 

Take some values of v and V = L n , and assume that the free energy satisfies 

X- 1 4 X J(T) = uf kB (T,v) + o(X- 1 ), 

as A —7- oo, where v := V/v is the volume per mole, and fk B (T,v) is a smooth 
function. In the limit A — > oo we have: 

df kB (T,v) p df kB (T,v) ^ 



&T ' dv 

where s — S/v, S is the entropy of the system, and p is the pressure in the system. 
This is precisely what one can see in the phenomenological thermodynamics. The 
equations (T5]) define a Lagrangian manifold L C 1R 4 (T, v, s, p) with respect to the 
symplectic structure 

uj := ds A dT — dv A dp. 

According to [H [TP] , the intensive thermodynamic coordinates (T,—p) can be 
perceived as "coordinates", and the extensive thermodynamic coordinates (s, v ) can 
be perceived as "momenta". More generally, one may perceive an abstract ther- 
modynamic system with d degrees of freedom £ = (£1,^2, •• • ,£d) as a Lagrangian 
submanifold A C M^ d with respect to symplectic structure u = J2i=i ^Vi A ^£i> 
77 = (771,772, •• ■ The Lagrangian manifolds are a central object in the semiclas- 
sical approximation of quantum mechanics [Tl] [12] and it is of interest to investigate 
the corresponding analogy between thermodynamics and mechanics in more detail. 
Recently this topic has attracted some additional attention and has received some 
new interesting developments in [31 [TS1 [HI HHJ 120] . 

In the present paper we prefer to use a little different coordinates to describe the 
thermodynamic Lagrangian manifolds. Denote e = Ejv the internal energy of the 



MICROLAGRANGIAN MANIFOLDS AND FLUCTUATIONS 



3 



system per mole. Then at temperature T = /3 1 we have: 

s = ^(f3f kB (f3,v)). 

The first law of thermodynamics tells us: 

ds(e, v) = fide + pdv, 

where p := p/T, and the entropy per mole is written as a function of e and v, 
s = s(e,v). The derivatives (3 = (ds/de) v and p = (ds/d(3) e define a Lagrangian 
manifold A thermo C M 4 (e, v , (3,p) with respect to the symplectic structure 

Uthermo ■= d(3 A de + dp A dv . 

Therefore, conceptually, the entropy (per mole) of a system is an action on A^ ermo . 
The total entropy S t hermo '■= vs is measured in the same units as the Boltzmann 
constant: 

[Sthermo] = [&b]- (3) 

If we have an abstract mechanical system with coordinates Q = (Qi,Q2, ■ ■ ■ ,Qn) 
and momenta P = (Pi, P 2 , . . . , P n ) and consider a Lagrangian manifold A mech in the 
phase space {R^ p ,u mec h), where 

n 

Umech '■— dPi A dQi, 

i=l 

then we have an action S mec h := J^PdQ, 7 is smooth curve on A mec h, which is 
measured in the same units as the Planck's constant: 

[Smech] = [h]. (4) 

In this context, the Heisenberg's uncertainty principle is totally similar to the Ein- 
stein's formula for the quasithermodynamic fluctuations: 

APiAQj ~ h, AEiA/3 m ~ k B , (5) 

where A denotes the standard deviation describing a fluctuation, i, j = 1,2, ... ,n, 
and I, m — 1, 2, and we put E\ = E, E 2 = V, fi\ — (3, (3 2 — P- Informally speaking, 
h is a quantum of mechanical action, and k B is a "quantum" of entropy 

The analogy between the Planck's constant h and the Boltzmann's constant k B 
expressed, in particular, by ([3]) , (j4]) , ([5]) , is the central motive of the present paper. 
The general scheme is as follows. First, we generalize the idea of quantization of phe- 
nomenological thermodynamic Lagrangian manifolds to the nonequilibrium setting. 
This defines a generic shape of the equations describing the transfer of fluctuations 
along a nonequilibrium evolution curve. Under some assumptions one can assemble 
these fluctuations into a "quasithermodynamic wavefunction" and extract the ther- 
modynamic Hamilton- Jacobi equation in the limit A — > 00. It turns out, that in 
nonequilibrium quasithermodynamics one can mimic Bell's inequalities and it is nat- 
ural to perceive the thermodynamic forces X = (Xi,X 2 , . . . ,X<i) as "coordinates", 
and the corresponding flows J = (J±, J 2 , ■ ■ ■ , Jd) as "momenta". We apply the sec- 
ond quantization to produce the "thermocorpuscles" introducing the creation and 
annihilation operators a ± (X, J) in a phase space point (X, J) G if^, the symplec- 
tic structure is := ^=1 dJi A dXi. Basically, the thermodynamic analogue of the 
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Wigner's equation for the single particle quasiprobability distribution is lifted to the 
symmetric Fock space J 7 *, 

T* : = C © L 2 (R™j) © (^(R^j))®**™" 2 © . . . , 

where <S> S ymm denotes the symmetric tensor product. After that, it is suggested to 
deform the second quantized equations by an interaction between the thermocorpus- 
cles and to mimic the basic constructions of statistical mechanics like the BBGKY 
chain and the collision integral. In particular, this yields a new interpretation of the 
nonequilibrium analogue Z^ X J(u;t) of the partition function in terms of the density 
of thermocorpuscles, 

Z { k X J(u;t) = J dX exp(-uX/k B ) J dJ (R[, a + (X , J)a~(X, J)R[), 

where u = (ui, u 2 , ■ ■ ■ , Ud) is a parameter varying in a neighbourhood of = 
(0,0, ...,0), R\ G J 7 * is a real vector depending on the rescaling parameter A 
and the time t, and (— , — ) denotes the inner product. 

2. Equilibrium quasithermodynamics 

Let us look at a phenomenological one-component thermodynamic system with 
two degrees of freedom, d = 2. To be more concrete, let the molar extensive coor- 
dinates be £i = e and £2 — v, where e := E/u, v := V/u, E is the internal energy, 
V is the volume of the system, v is the number of moles. We have a phase space 
R 4 (£i) £2, V2, 772) with a symplectic structure oj = dr/iAd^i + d^Ad^, where rji — j3 is 
the inverse absolute temperature of the system, and r] 2 = p := (3p, p is the pressure 
in the system. The system is described as a two-dimensional Lagrangian manifold 
A in this phase space defined by equations 

ds\ „ /ds \ 
de)v ^ \dv)e 

where s — S/v is the entropy S of the state of the system taken per mole. The 
function s = s(e, v) is a generating function of A with respect to the focal coordinates 
(s,v). If we assume, that A admits another choice of focal coordinates, say, (/3,p), 
then the corresponding generating function Jl = jl(/3,p) is linked to s = s(e,v) via 
the Legendre transform, 

K^P)\ k = [s(e,v) -pe-pv]\ A . 

The physical meaning of this function is just Jl = —(3fi, where fi is the chemi- 
cal potential of the system. It is natural to consider along with A a Lagrangian 
manifold A + of dimension d + 1 = 3 in the phase space M. 6 (x,y) with coordinates 
x = (x 1 ,x 2 ,x 3 ), xi = E, x 2 = V, x 3 = v, and "momenta" y = (y u y 2 , y 3 ), y 1 = (3, 
Vi = Pj U3 = A*- The symplectic structure u + is given by u + := Y^=i dy% A dxi, and 
the manifold A + can be described as 

A + := {(x,y) \x 3 ^ 0, (x 1 /x 3 ,x 2 /x 3 ,y 1 ,y 2 ) e A, y 3 = ji(x 1 /x 3 ,x 2 /x 3 )}. 

One can observe, that A + admits a projectivization with respect to the extensive 
coordinates x = (xi,x 2 ,x 3 ) = (E,V,u), which can be taken as global focal coordi- 
nates on A + . At the same time, all the points of A + are singular with respect to the 
focal coordinate plane M 3 (y), y = (yi, 2/2, 2/3) = (/3,p,jj)- Our physical system can 
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be identified with this A + , and one should assume (for the physical reasons), that 
A + is connected and simply connected [8j [H] . 

It can happen, that A + admits other focal charts with a choice of coordinates 
different from x = {x\,x 2 ,x 3 ). Recall the notation: 

(aci,a5 2 ,a? 3 ) = (E,V,u), (2/1,2/2,2/3) = (P,+Pp,~Pfi). 

In a nonsingular chart U with coordinates (x\, X2, x 3 ) = (E,V,v), the generating 
function S(xi,x 2 , x 3 ) is the entropy of the system in a point a G U C A + correspond- 
ing to (xx, x 2 ,x 3 ), (dS - Y^=iyidxi)\A+ = 0, and S(E,V,v) = us(e,v), e = E/v, 
v = V/v. In a singular chart U\ with coordinates (j/i, x 2 , x 3 ) = (ft, V, v) the generat- 
ing function Si (2/1, x 2 , x 3 ) can be interpreted (up to an additive constant) as —/3F, 
where F is the free energy of the system in a point a G U\ C A + corresponding to 
(yi,x 2 ,x 3 ). In a singular chart Ui )3 with coordinates (yi,x 2 ,y 3 ) = (/3, V,— PfM), the 
generating function 61,3(2/1, x 2 , 2/3) can be interpreted (up to an additive constant) 
as — /3f2, where Q is the thermodynamic potential with respect to (T, V, /1) taken in 
a point a G C/i,3 C A + corresponding to (yi,x 2 ,y 3 ). 

Let us analyse the case of the to charts U\ and U\ i3 in more detail. Assume 
U\ fl £/i ; 3 7^ and fix a G f/i fl C/i )3 . In the ambient space M?(x,y) D A + this point 
acquires coordinates x(a) = (xi(a),x 2 (a),x 3 (a)) and y(a) = (2/1(0), y 2 (a), y 3 (a)). 
The step from phenomenological thermodynamics to statistical thermodynamics cor- 
responding to the first chart consists in the following. One considers a family of 
physical systems parametrized by A > (the rescaling parameter), A — > 00. Ev- 
ery system is placed in a thermostat and the walls of the system are fixed and are 
impenetrable for the particles. The inverse temperature P\(a) = x±(a) is the same 
for each A, the number of particles N\(a) = \x 3 (a)R/kB, together with the volume 
of the system V\(a) = \x 2 (a), grow linearly as A — > 00. For simplicity, one may 
have in mind a system of particles moving on a n- dimensional torus of radius L 
described by a Hamiltonian H^\g), where g is a collection of parameters describing 
the interaction between the particles and the geometrical shape of the volume (the 
external potential). One substitutes L = (V\(a)) 1 ^ n and A^ = N\(a), assuming that 
the parameters g = g\(a) may be adjusted as well. The free energy F^ N (T; g) at 
absolute temperature T is given by F^ N (T;g) = —kBT\nZ^ N (T~ l ;g), where 

Z^ N {P;g) := Trexp ( - -l/3i^%)), 

is the partition function of the canonical Gibbs distribution at absolute temperature 
T = fi~ l (we assume that the corresponding trace is finite for the required values of 
(3, L, N, and g). The link with phenomenological thermodynamics is as follows: 

5i( yi (a),x 3 (a),.T 3 (a)) = hm [-A" 1 /^) ^n^^^ 9x(a))+C(X)], (6) 

where it is assumed that one can find functions g\(a) and C(A) so that the limit on 
the right-hand side exists. Denote S[ X \a) the expression in the square brackets on 
the right-hand side of 

Now if we take the other chart Ui^ with coordinates (y\,x 2 ,y 3 ) = ((3, V, — f3fi), 
then we need to consider another family of rescaled systems of another sort: the 
walls of the system must admit a penetration of particles, i.e. the border of the 
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system is described purely geometrically. One still has a rescaling parameter A > 0, 
the volume of the rescaled system V\(a) = \x 2 (a), but now one should look at the 
partition function ^j^(/3,p;g) of the grand canonical ensemble, 

(l L J ((3, ng):= Tr exp ( - ±-0 [fiM (g) - ^] ) , 

(it is assumed that the trace is finite for the required values of (3, /i, and g), where (3 
and fj, are the inverse absolute temperature and the chemical potential determined 
by the environment of the system, respectively, l-L^ L \g) is the second quantized 
Hamiltonian of the multiparticle system on a n-dimensional torus of radius L, g is 
a collection of parameters describing the interaction between the particles and the 
external potential, is the second quantized operator of the number of particles 
in the system j2]. Denote flj^ (T, fi; g) := — fcsTlnQ^(r _1 , p; g) (the potential Q). 
One needs to substitute in place of (3, p, and g some functions ft\(cx), p\(a), and 
g\(a), respectively, in such a way, that the average energy and the average number 
of particles in the system grow linearly with A — > oo. It is convenient to change the 
variables yx := /3, y 2 := — ftp, and to consider Z^(yx, y 2 ; g) := Cka (ft > I 1 '-, 9) ■ Denote 
(ft\(a; g), — ft\(az] g)fi\(a; g)) the solution of the system of equations 

(^ < "' ,1/ " ) (v..^;s))"(-^|;)^ <<,))I '" , fe.v 2 ; S ) = Ax 1 ( a ), 

with respect to (yx, y 2 ), assuming it exists and is unique for the required values of g. 
Recall, that R denotes the universal gas constant, Xi(a) corresponds to the internal 
energy, x 2 (a) corresponds to the volume, and x 3 (a) corresponds to the number of 
moles of the chemical substance in the phenomenological system. The link with the 
phenomenological thermodynamics is as follows: 

Si,3(yi(a),x 2 (a),y 3 (a)) = 



lim 

A— >oo 



A-%(a^)4 ( f 2(a)rl) (/3 A (a^)-\ yUA (a^)^)+(7(A) , (7) 



9=g\{a) 



for some functions g\(a) and C(A) ensuring the existence of the limit on the right- 
hand side. Denote S'j^a) the expression in the square brackets on the right-hand 
side in (j7|). 

The two functions Sx(yx, x 2 , x 3 ) and Sx^(y\,x 2 ,y^) are related via the Legendre 
transform in the third argument. If we perceive them as functions on the thermo- 
dynamic phase space IR 6 (a;, y), then we have 

Sx,3(yi,x 2 ,y 2 )\ A+ (a) = [Sx(yx, x 2 , x 3 ) - y 3 x 3 ] \ A+ (a), 

for any a G Ux D U\$. Since A + admits projectivization with respect to the ex- 
tensive coordinates x = (xx,x 2 ,x 3 ), the generating functions in the focal charts are 
homogeneous with respect to these coordinates: 

Sx(yi, px 2 , px 3 ) = pSx(yx,x 2 ,x 3 ), S 1>3 (yx, px 2 , y 3 ) = pS 1<3 (yx,x 2 ,y 3 ), 

for any p > 0. 
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Consider now an abstract phenomenological thermodynamic system with exten- 
sive coordinates x = (xi,x 2 , . . . , Xd+i) and intensive coordinates y = (y%, y 2 , ■ ■ ■ , Vd+i)- 
The corresponding thermodynamic phase space is defined as M 2< - d+1 ^(x, y) equipped 
with canonical symplectic structure cu + = J2i=i dyi A dxi. The equilibrium states of 
the system are identified with a Lagrangian manifold A + C R 2 ^ +1 -*(x, y), which is 
connected, simply connected, contained in the region Xi > 0, % = 1, 2, . . . , d + 1, 
and admits a projectivization with respect to x: if (x(a),y(a)) are the phase 
space coordinates corresponding to a point a G A + , then for any p > 0, the co- 
ordinates (px(a),y(a)) correspond to another point in A + . Let / be a subset of 
[d+1] : = {l,2,...,d + l}. Denote 

(x, y)/ := ({x j } je[d+1] \ I , {y^iei) , (8) 

where the indices atop denote the order of appearance of the variables in the list read 
from left to the right, for example, if d = 2, I = {1,3}, then (x,y)i = (y±, x 2 , 2/3). 
In particular, for the emptyset / = 0, we have (x,y)q) = (xi,x 2 , ■ ■ ■ ,Xd+i), and for 
I — [d + 1] we have (x, y)[d+i] = (2/1, 2/2, • • • , Ud+i)- A focal chart of type / on A + is 
an open set U C A + together with a fixed choice of local coordinates of the shape 
(x,y)j. For every focal chart of type I over U, there exists a function Suj((x,y)i), 
such that U is described by the equations: 

_ dSu,i((x,y)i) dSujijx^y)!) 
dxi dyj 

where i G /, j G [d+ Observe, that for any p > 0, 

Su,i((px,y)i) = pSuj((x,y)j), 

whenever both left- and right-hand sides are defined. Observe also, that the ther- 
modynamic Lagrangian manifold A + never admits a focal chart of type + i.e. 
at least one of the local coordinates in a given focal chart must be extensive. In 
case U is connected, the function Sjj,i is defined up to an additive constant, and for 
a fixed pair of focal charts (Ui,I) and (U 2 , J), such that XJ\ fl U 2 ^ 0, it is always 
possible to adjust these constants in such a way that 

Su 2 ,j((x,y),j) = S Ult i{[x,y)i) - ^ yjXj + ^ 

jeJ\i iei\J 

for every (x,y) corresponding to a point a G U% fl U 2 . This implies, that Su u i and 
Su 2 ,j are linked via the Legendre transform (denote it C U j 2 P Ul ), Su 2 ,j = ^ l j 2 p Ul {Su 1 j), 
for every (x,y) corresponding to a point a G Ui fl U 2 . Observe that £f j = id, for 
a focal chart (U, I), where U is connected. If U admits focal coordinates of another 
type J, then we have CJ{ J oCj J J1 = id. Furthermore, if U admits focal charts of types 
/, J, and K, then it is straightforward to check the cocyclicity condition: 

L\ K o C U KJ o = id. (9) 

Denote ttj : A + — > M. d+l the canonical projection a H- (x(a),y(a))i, for I C [d + 1]. 

Return now to the example discussed above, d = 3. We have defined the func- 
tions S[ (a) and S[ X ^(a) (see the definition right after the equations (jSJ) and ([7]), 
respectively), such that Si(yi(a),x 2 (a),x^(a)) = lim^oo s[ x \a), for a G Ui, and 
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Si (3/1(0;), x 2 (at) } 3/3(0:)) = limA-^oo £13 (a), for a G ?7 1)3 . What happens in the higher 
order corrections with respect to A -1 — > 0? In quasithermodynamics one assumes 
that these functions admit the following asymptotic expansions: 

00 00 
S[ X \a) ~ XJA-W), S$J(a) ~ X)A^g(a), 

Z=0 i=0 

where <£>i (a) and <£>i 3(0) are smooth functions, Z = 0, 1, 2, This can be perceived 

as a condition on a choice of the function g\(a) present in the definitions. It is a 
little more convenient to work in terms of coordinates, introducing the functions 
,x 2 ,x 3 ) and $13(2/1, #2, 2/3)5 defined as follows: §f\yi(a),x 2 (a),x 3 (a)) = 
if? {a), a G C/i, and { Vl (o), x 2 {a), 2/3(0')) = rf^o'), a' G U 1>3 , I G Z >0 . We 

know that $S° } (?/i, ar 2 , x 3 ) = S 1 (yi,x 2 ,x 3 ) and $!°J(2/i, x 2 , y 3 ) = S lt3 (y u x 2 ,y 3 ), so 
the link for the leading coefficients is the Legendre transform. What is the link 
between the higher order coefficients corresponding to I ^ 1? 

It is more convenient to describe the link mentioned in an abstract setting. Con- 
sider a Lagrangian manifold A + C M 2 ( d+1 ) (x, y) for an abstract phenomenological 
thermodynamic system as above. We assume that in every focal chart (U, I) of type 
/ C + we are given a function S^\(a), a G U, depending on a parameter A > 0, 
such that there is an asymptotic expansion 

00 

4>)^E A ^(«), ( 10 ) 
1=0 

as A —7- 00, where (p^j(a) are smooth functions, I G Z^q- Using the notation (jSJ), 

define also the functions $^j((x, y)/) via $^ 7 ((a;(a), y(a))i) = </?^j(o), a G U, for 
I = 0,1,2,.... Take another focal chart (W, J) on A+ such that W n U ^ 0. In 
the leading order / = 0, in the points corresponding to o G W D £/, we have the 
Legendre transform: 

The formulae for the other coefficients mimic basically to the formulae of the sta- 
tionary phase method. We need to assume that the Hessian 

detHess^ ) / ((x,2/)/)^0, 

for every (x,y)i G tti(U), for every focal chart (U,I). 

Let us first introduce some notation for the differential operators associated with 
this method. If we have a function f((x,y)i) G C^injiW fl [/)), then we can look 

at an integral Qwr\U J A^un /](( x > v)j)i £ nj(W fl U) depending on a small 

parameter e > 0, 

Q!^XVK(^) : = (2n£) \Ai\ + ii\J\ I ( II <*»*)( II dy l )f((x,y) I )x 

^ ' J j£J\i ieJ\J 

x exp |i $^((a;,2/)/) - ^ + 2/i%] }, 



ie/\J jeJ\/ 



MICROLAGRANGIAN MANIFOLDS AND FLUCTUATIONS 



9 



where | • | denotes the cardinality of a set. This integral admits an asymptotic 
expansion as e — > (the stationary phase method, see [5]): 

oo 

x exp {^flj") E(- i£ ) n ^i[<'/](^(^^)^))' 

n=0 

where the map Tfp v '■ kj{W fl U) — > tti(W fl U) is the glueing map between the fo- 
cal charts (U, I) and (W, J) on A + defined by r]^ 7 nf/ ((x(a;), y(a)) j) = (x(a),y(a))i, 
for a E W C\ U, /ijj- nC/ € Z/4Z is a constant (related to the Maslov index), and 
^wnu jA^uii ~ ] are nnear partial differential operators, such that, for each n, 
Vwr\U jA^uii f](( x > v)i) depends only on a finite number M n of partial derivatives 
of <&u\ in the point ((x, y)i) G tt i (W fl £/). The number M n grows as n — > oo, and in 

the leading order n = the operator V^ nU j j[^u\, — ] is just a multiplication over 
a smooth function, which does not vanish in no point of nj(W D ?7). Fix the choice 
of i n such a way that this function is positive (this is always possible). 

Observe, that if we take instead of f((x, y)j) a function f( e '((x, y)j) depending on 
the small parameter e, such that there is an asymptotic expansion f^ e '({x,y)i) ~ 
Y2=,{-^) m fm{{x,y)i), where f m G C^i(W fl U)), m = 0, 1,2, . . . , then, using 
the linearity of the operators V^ nUJI [^ l7 — ], n = 0, 1, 2, . . . , we obtain: 

ex? (- ^Z(^yh) - j/*Sr n,7 )^w[^i,/ (e) ]((*.y) J ) * 

oo n 

-E(- i£ ) n Ett[<>/ m ](^T((^)^)). (12) 

n=0 m=0 

for (x,y)j G ttj(W fl Z7). Now, let us assume that p E > can be written in the form 
= exp{g^}, where g^ is a convergent power series g^ = Ylm=o(~ i £ ) m 9m/m\, 
where g m G C^(jti(W fl (/)), m = 0,1,2,.... Assume also that the left-hand 
side of (fl2~|) (denote it at this moment F^ £ \(x,y)j)) can also be represented as 
F( £ )((x, y)j) = exp {G( e )(r/j nC/ ((x, y)j))}, where G^ e \{x,y)i) is a convergent power 

series C7( e ) = E"=o(- i£ ) mG V m! . with G ™ e Cg°(7f/(W n 17)), m = 0, 1, 2, . . . (at 
least this does not contradict a possibility of having an asymptotic expansion like 
on the right-hand side of (fl2l ). The link between the collections of the coefficients 
{G m }m=o an d {dm}m=o defines some operators, that we are going to use to describe 
the link between the asymptotic power series (ITU]) corresponding to different (U,I). 

Take a pair of smooth functions a(£) and &(£), £ varies over R, a(£) = exp(6(£)). 
Look at the Taylor's expansions a(£) = X^m=o am ^ m / m ' anc ^ = S^Lo ^n£"/ n '- 
For the leading coefficients, a = exp(6 )- If m ^ 1, then 
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m ^ 

b \h\b\h\---b\h\- 

P =1 h,I 2 ,-,I P C[p], 



h,I 2 ,...,I p ^0, 
\h\ + \h\+---+\I P \=m 



Therefore, a m = A m (b ,b 1 ,b 2 , . . . , b m ), where 



A m (b ,b 1 ,b 2 , • • • ,b m ) := e b ° S~] — — - Tb mi b m2 ...b m 

p=l J mi,m2,...,m p =l, " ^ 

mi+rrt2H hwip=m 

(13) 

for every m = 1, 2, To invert these formulae, observe that b = lna , a > 0. 

For m > 1, we have 



6 - = {(IP^ML = {(^^^KO/ao - 1]'} 



9 (-1) 



9=1 



5=0 



= — q- y~) a \ji\ a \.h\ ■ ■ ■ a \j q \- 

9=1 ^ U Ji,J 2 ,...,J 9 C[ g ], 

\Ji\ + \J2\ + ---+\J q \=m 

Therefore b m = B m (a , ai, . . . ,a m ), where 

v ^ (— l) 9 ^ m ' 

-t'ml^O) ^1) • • • j Q>m) • / , a / , ! i r^mi ^m2 • • • ®"m a i 

qak *— ' mi!m 2 !...m ! 

9=1 Hll|TO2,...,»T(g=l, 
m\+m2-\ hmij=m 

(14) 

for every m = 1, 2, . . . . Extending naturally the notation as 

:= exp(6 ), B (a ) : = ln(a ), (15) 

we arrive at a n = v4 n (6 , &i, • • • , b n ) and b n = B n (a , oi, . . . , a n ) for all n = 0, 1, 2, 

We can now describe how the expansions ffiUj) corresponding to different focal 
charts (U,I) on A + must be linked. Take a pair of focal charts (U,I) and (W, J), 
W fl Z7 7^ 0, and look at the expansions for S'^j(a) and S 1 ^ j(a) in terms of the local 
coordinates: 

oo oo 

where (x, y) corresponds to a point a varying over WHU, and S^foy)/) := 
5 'S(^7 1 (( X )2/)/)) for (^1/)/ e ^/(f 7 ), and y)j) := S&C^Cfo y) j)), for 



(x,y)j G 7Tj(W). To compute the collection of coefficients j(( a; ) 2/) j)}j=o- 

{x,y)j G 7rj(W fl Z7), from a collection of coefficients y)/)}^o> one nee ds 

to do the following: 

• Compute &wj(( x iV)j) as a Legendre transform (JIT]) of $^ J ((x, y)/). 

• Set 6 m = y)i), m = 0,1,2,..., and compute the corresponding 
a n = A n (bo, b n ), n = 0,1,2, ... , where A n is defined in ffl3|) and f fT5|) . 



I oo 
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• Compute the coefficients of the asymptotic expansion on the right-hand side 
of (H2D, replacing {f m }™=o with { a m}m=o> a ™ = a m ((x,y)i)- 

n 

for n = 0,1, 2,.... 

• Compute b' n := B t ,(a' 11 4,...,oJ l ), n = 0,1,2,..., where B n is defined in 
f }T4|) and ffT5j) . This collection of functions b' n = b' n ((x,y)j) is precisely the 
required collection of the higher order coefficients, 

where n = 0, 1, 2, . . . , and (x, y)j corresponds to a point in W fl U. 

It can be of interest to assemble the higher order coefficients ?/)/), / = 

1,2,3,..., into a formal power series 

oo 

where e is a formal variable, $$\((x, y)i) G C7°°(7rj(J7))[[e]]. Note, that the coef- 
ficient $^j((a;, y)/) is noi involved in this definition. If U admits another type of 
focal coordinates J C [d + 1], then the described link between the coefficients in 
different focal charts yields a nonlinear map 

We have M^j = id, and, in case U admits the focal coordinates of types J, J, and 
K, one can check the cocyclicity condition: 

^"Af.VAfJ^i (16) 

Note, that about U we assume that the Hess's matrices HessSj and HessSj are 
non-degenerate. Intuitively, the map Mj j "extends" the Legendre transform Cjj, 
and the cocycle condition (fl6|) corresponds to the cocycle condition Qj. 

Where do the functions S^\(a) come from? Consider again a phenomenologi- 
cal thermodynamic system with extensive coordinates (xi,X2,xs) = (E,V,u) (in- 
ternal energy, volume, number of moles), and intensive coordinates (2/1,2/2,1/3) = 
(/3, (3p, —fill) (/3 is the inverse absolute temperature, p is pressure, \i is chemical po- 
tential). If we know the entropy of the system as a function of extensive coordinates, 
S = S(E, V, u), then the first law of thermodynamics can be expressed as follows: 

3 

dS = (3dE + (f3p)dV + {-^)dv = Vi dx h ( 17 ) 

i=l 

on the Lagrangian manifold A + . If / = {1} and U C A + is fixed, then the limit 
S\(a) := limx^oo S^](a) corresponds to the Legendre transform of S in the first 
argument, dS\ = (—Ed(3 + ((3p)dV + (— (5jj)dv)\u, i.e. Si = —/3F, where F is the 
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phenomenological free energy described as a function on U C A + . If one takes the 
canonical Gibbs distribution 

u ^ (ft9) "4^ exp (-^H- (18) 

then one can see that the averages (H^ ) and ((H^' — (H^)) 2 ) over w[ L J N (P;g) 
correspond to the derivatives {—kBd/dfi) 1 and (— &;£<9/<9/3) 2 of hxZ^ N (/3; g). Sub- 
stituting the values (5 = (3\(a), L = (V^ct)) 1 /™, N = N\(a) = u x (a)R/k B , and 
g = g\(a), corresponding to the rescaled system (the rescaling coefficient A), one 
extends this fact as follows: the quantities 



(/3,V,v)=(f3 x (a),V x (a),v x (a)) 



define the cumulants of the fluctuations of energy, n = 1,2,3,..., where 7?i : 
a i — Y ([3 (a), V (a), is (a)) is the canonical projection. Observe, that C[ X \a) = 0(A), 
A — > oo, which is the expectation value of internal energy, and C { 2 X \a) = 0(A), the 
variance of the corresponding fluctuations, so the standard deviation is 0(y/X), as 
it should be in quasithermodynamics. The inverse absolute temperature (3 is fixed 
by the thermostat, and the canonically conjugate extensive quantity (the internal 
energy) fluctuates. This interpretation is naturally dualized: if an extensive quantity 
is fixed (for example, the volume V), then the canonically conjugate intensive quan- 
tity fluctuates. From the first law of thermodynamics ( jl7JL we can see, for example, 
that the derivatives over V should define the cumulants Cm\o) of the fluctuations 
of (3p, where p is the pressure in the system. More precisely, the derivative d/d(3 
should be replaced with X^d/dV, since V\{ot) = 0(A) (extensive variable), and 
/3\(a) = 0(1) (intensive variable), A — > oo, and we should take into account the 
minus sign corresponding to the Legendre transform: 

v ov 7 Kb {py,u)={p x { a )y x { a ),v x ( a )) 

for m — 1,2,3, In particular, the first two cumulants for the fluctuations of 

(3p (the expectation value and the variance) satisfy C[ X \a) = 0(1) and C 2 X \a) = 
0(A _1 ), A — > oo, so the standard deviation corresponding to f3p is 0(l/yfX), as it 
should be in quasithermodynamics. 

For an abstract phenomenological thermodynamic system A + C IR 2 ( d+1 ) (x, y) the 

interpretation of the functions S^\(a) for every focal chart (U, I) is as follows. We 
know that [d + 1]\J cannot in thermodynamics be empty (by construction). 
• The first law of thermodynamics is expressed over (U, I) as follows: 



d fuj = ( ~ Xidyi + y i dx > 



iei je[d+i]\i 



where (p^\(a) := lim^oo S^](a), a e U C A + , x — (xi, x 2 , ■ ■ ■ , Xd+i) are the 
extensive coordinates, and y = (yi, y 2 , ■ ■ ■ , Vd+i) are the canonically conju- 
gate intensive coordinates. A model example is d — 2, (xi, x%, X3) = (E, V, v) 
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(internal energy, volume, number of moles), (1/1,1/2,1/3) = (P, Pp, — P/i>) (in- 
verse absolute temperature P, pressure p, chemical potential /i) 

• The units of measurement of Syl(a) are same as the units of measurement 

of the Boltzmann constant ks, [Sui(a)] = [ks]- The function (p^\(a) corre- 
sponding to the case / = is the entropy of the phenomenological system in 
the state a G U C A + . 

• The following collection of derivatives is interpreted as the cumulants de- 
scribing the fluctuations of the vector (x, y)[d+i]\i over a state a G U: 

r (U,I,\) ( x| _ 

°m K a )\ a= ^ l((x , y)l) — 

rK-^rx n (^n ^-'» . 

iei y ie[d+i]\/ J 

where M = (mi, rri2, ■ ■ ■ , Wd+i) G Z^, 1 is an integer multi- index, M 7^ 
(0, 0, . . . , 0), and ttu,i '■ U 3 a (->■ (x(a), y(a))/ denotes the canonical projec- 
tion, (x,y)j G 7fi(U). The cumulants corresponding to the same 
but computed over two different focal charts (U, 7)and (W, J), (i.e. different 
families of rescaled systems) coincide asymptotically in the intersection: 

C^' X) (a) = C^ x \a) + 0(X-n, 

as A —¥ 00, where a G U D H 7 . 

The cocycle conditions ()9]) basically say that one can glue a Lagrangian manifold 
A + . Informally speaking, phenomenological thermodynamics IS a Lagrangian man- 
ifold A + (connected, simply connected and admitting projectivization with respect 
to extensive coordinates). The statistical thermodynamics IS the Gibbs distribution 
corresponding to a focal chart (U, I) on A + , / C [d+1], I ^ [d + 1], for example, the 
canonical distribution ( ITS]) if J = {i } is a singleton and i corresponds to the inverse 
absolute temperature. The type I of the focal chart describes the way we "extract" 
the system from the outside world (fixed adiabatic walls, fixed heat-conducting walls, 
a non-fixed piston, walls penetrable to certain sorts of particles, etc.) This way of 
extraction does matter if we consider the Gibbs distribution along the families of 
rescaled systems by a parameter A — > 00. If one is ready to sacrifice the precision 
modulo 0(X~°°) describing the fluctuations of the intensive and extensive quantities, 
then one arrives at "quasithermodynamics" . One may say that quasithermodynam- 
ics IS a global section of the sheaf on A + corresponding to the cocycle conditions 
(|T6]) (the formal variable e corresponds to the small parameter A -1 , but should not 
be confused with it). This sheaf corresponds to the sheaf of V-objects in [12], and the 
analogy with mechanics is as follows [8]: phenomenological thermodynamics corre- 
sponds to classical mechanics, statistical thermodynamics corresponds to quantum 
mechanics, and quasithermodynamics corresponds to semiclassical mechanics. 

3. NONEQUILIBRIUM QUASITHERMODYNAMICS 

We are now interested in extending the quasithermodynamics picture to a nonequi- 
librium setting. One may use the concept of relevant ensembles for this purpose |23j. 
Like in the previous section, to keep the story more simple, consider a multi-particle 
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quantum mechanical system on a n-dimensional torus of radius L. Denote T the 
Fock space of the system, and % the Hamiltonian of the system (a self-adjoint oper- 
ator on J 7 ). Fix a collection of observables S[ L \g) , 8^\g) , . . . ,S d L \g) (self-adjoint 
operators on the Fock space of the system), depending on a collection of parameters 
g = (<7i, #2, ■ ■ ■ , g m )- One considers the Gibbs's statistical operator 

where zjj. (/?; (?) := Tr( — Y^t=i fii^i (sOAs) (we assume that the trace exists), 
is the Boltzmann constant, and (3 = (/3j, /3 2 , . . . , /3^). Set 

The general idea is that the evolution of the system is well approximated by a family 
of operators {w^ \f3 f , <?*) }., where t is time, and L*, /?*, are certain functions 
(we let L and g change with time for more generality). One says that the state 

Tr{{yj^ <?*) — } defines a relevant ensemble at the moment of time £. 

Let us look at the collection of all possible relevant ensembles, i.e. all possible sta- 
tistical operators w^J(/3;g). To mimic the setup of the equilibrium thermodynam- 
ics, consider a phase space with symplectic structure u + := Ylt=i dyi A dxi, 
x = (xi, x 2 , ■ ■ ■ , a?d+i), y = (yi, V2, ■ ■ ■ 1 The first d coordinates correspond 
to Sj; L \g), i ^ d, and the last coordinate Xd+i corresponds to the volume V = L n . 
The first d coordinates correspond to i ^ d, and the last coordinate yd+i 
can be perceived as a non-equilibrium analogue of the "pressure over temperature" 
variable in the equilibrium thermodynamics. The coordinates x are termed the ex- 
tensive coordinates, and the coordinates y are termed the intensive coordinates. 
Assume that there exists a Lagrangian manifold A + C Mx[y +1 \ which is connected, 
simply connected, admits projectivization with respect to the extensive coordinates 
(i.e. the analogue of the Lagrangian manifold in phenomenological equilibrium ther- 
modynamics). Assume, for simplicity, that A + admits (yi, y 2 , ■ ■ ■ , yd, Xd+i) as global 
coordinates (i.e. is covered by a single focal chart of type [d] C [d+ 1]). Let A — > oo 
be a large parameter (termed the rescaling parameter) , and assume that there exist 
functions C(A) and 0\(a), g\(a), where a varies over A + , such that 

OO 

X~ 1 ^^ Xd+1 (a))1/n) (/? A (a) ; g\(ai)) - C(X) ~ ]T A^>), (20) 

1=0 

where <f>^ B (oc) are smooth functions, I = 0, 1, 2, ... . This expansion should be per- 
ceived in analogy with (flUj) . and V\(a) = \xd+i(a) corresponds to the expanding 
volume as A — > oo. 

Like in the equilibrium case, we should draw a distinction between phenomeno- 
logical, statistical, and "quasi" nonequilibrium thermodynamics. The general idea is 
that the phenomenological nonequilibrium thermodynamics is described as a curve 
7 = on the Lagrangian manifold A + , where t is time. If this is a curve of 

relaxation to equilibrium, then some points of A + should correspond to the equi- 
librium states. As a model example, one may think of a system consisting of two 
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subsystems which interchange heat with each other, but have different temperatures 
(the volumes of the subsystems are assumed to be fixed). The internal energies E\ 
and E 2 of these two subsystems (indexed by 1 and 2) will change in such a way 
that Ei + E 2 is fixed, until the inverse temperatures j5\ and (3 2 become equal. More 
generally, if, for example, the total volume of the system V(a) = x d +i(oi) is fixed 
along the evolution 7 = {a*} t , Vi^a 1 ) = const, then under some assumptions, like 
"the smallness of time memory effects" [23], one can derive the following system of 
equations for Ei(t) := ^(a*), i — 1, 2, . . . , d: 



dEi{t) 



JiW+^LijWPjQ), (21) 
3=1 

where (3i(t) = ?/«(«*), i G [d], Lij(a) are the Onsager coefficients corresponding to 
a G A + (in practice, they are often approximated as constants) , and Jj(a*), i G [d], 
are the flows induced by {Pj(t)}je[d\- These flows can be perceived, for each i G [d], 
as a limit 

Ji(a):=]hn (\- 1 TrfwP(B]g)j} L \g)} ), (22) 

where J^ L \g) are the operators of flows, 

jl L \9)--=\\H( L \g)A L \9)\ (23) 

the square bracket denotes a commutator, and h is the Planck's constant. 

Denote the left-hand side of ( 1201) as sj£j(a), and let S^(yi,y 2 , ■ ■ ■ , yd, %d+i) be 
this function expressed in terms of the focal coordinates on A + of type [d] C + 

S { k X g(yi(a),y2(a),- ■ -,yd(a),x d+1 (a)) = S ( k ^(a) } 

for a G A + . In analogy with the previous section, one interprets the derivatives 

<*>(■> - ((n (-<n K^r ^r^ 1 ' ) L<»>- 

(24) 

where a G A+ C Af = (mi, ...,m d , m d+1 ) G Zg 1 , M ^ (0, 0, . . . , 0), as the 

cumulants of fluctuations of (xi, . . . , x d , yd+i)- The Lagrangian manifold A + is re- 
covered from the leading term . . . 2^+1) := lin^oo Sj^iyi, ■ ■ ■ , y d , x d +i) 
of the asymptotic expansion in A -1 as a system of equations: 



where i = 1, 2, . . . , d. 

It follows that in every point a 1 of the curve of the evolution 7 = {a t } t , we have 
a collection of cumulants {C^ (ck*)}m describing the fluctuations of the measured 

values of the quantities associated with s[ L \g) , S 2 L \g) , . . . , S d L \g) and p, where 
p is the intensive quantity dual to the volume (the nonequilibrium "pressure over 
temperature"). This is, of course, an approximation, depending, essentially, on a 

successful choice of the basis variables £- L \g), i G [d], for the family of relevant 
ensembles. These collections of cumulants stem from the same Lagrangian manifold 
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A + , but, in fact, one can be more general, and consider germs of different Lagrangian 
manifolds "attached" to each point a 1 £ 7 C A + C R^ +1 ^. To see this, one 
needs to consider the generalized Fokker-Planck equation [24J. The idea of the 
derivation of this equation is totally similar to the one for ( J2TT) . except that instead 
of a finite collection of basis quantities {£} (g)}ie[d]i one considers a continuous 

collection {L n 8 £ {a x - £^'{9)) ■ ■ ■ S e {a d - S^ig))} (ai,...a d )eM d ) where 5 £ (-) is a smooth 
approximation of the Dirac's delta, 5 e (-) —> 5(-) as e — > (in the weak sense), 
J da5 £ (a) = 1, a = (01, 02, ■ ■ ■ , a<z)- Of course, a special care needs to be taken about 
treating a delta function of an operator, and also there is a technical problem related 
to the noncommutativity of the operators £j, % £ [d]. Let us assume, for simplicity, 
that S[ L \g), 82 (g), ■ ■ ■ ,Sj i L \g) mutually commute (the noncommutativity can be 
handled as well, see, for example, [16J). Then, for every a = (ai, a 2 , ■ ■ ■ , a^), we can 
just write 

d 

^ L \a ] g):=L n \[8 £ {a i -£l L \g)), 

i=l 

without paying attention to the order of operators. Assume that the corresponding 
operators of flows j[ L \g) , J2 L \g) , . . . , J^ L \g) defined in (1231) mutually commute 
as well. 

In place of {S- L \g)} iG \d], we have a collection {T £ (a; g)} a m d > but with an addi- 



Ue[d] 

tional condition 



J da^ L \a;g) = L r ' 



One must also point out, that all Fe L \a',g), a £ M. d , have the same units of 
measurement. In this sense, {T £ L \a\ g)} ae u d corresponds to {£} L \g)}f =0 , where 
£y L \g) := L n — X/i=i ^ (flO 5 if we choose all £- L \g), i £ [d], having the same units 
of measurement as the volume L n , so that we have a right to add them. The Gibbs's 
statistical operator ( !T9|) can be expressed as follows: 

s)= exp / _ 1 m^^ (9) + ±{ Pt+ y 



Denote 



Wg{m 9) ■= exp ( - i £ t/^Q/)) , (25) 

where 77 = (770, 771, . . . , 77^) is a collection of parameters. Observe that the trace 
of this operator is equal to one if we substitute 770 = L~ n &^ (f3; g), and 77$ = 
fii + L~ n $j^(/3; g), i £ [d\. One can also check, that 

"' d 



(E^) lnT ^>;s)) = ^ 

i=0 /l 



for any 77 = (7/0,7/1, . . . ,7/ rf ). Consider a phase space with coordinates with 
coordinates 7/ = (770, 771, . . . , 77^), and £ = (£ , £1, • • • , £d)> an d symplectic structure 
a) := J2i=odVi A The map r : -)■ (£,77) ^ (2,7/), defined by 
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X 



Ci, Vi = Vi - Vo, for ^ % d, and x d+1 



symplectic. Therefore, in place of the Lagrangian manifold A + C M.x^ +1 ^ 



£o + Ei=i6, Vd+i = Vo, is 
i 

one can 

work in terms of the manifold r _1 (A + ) C K.^ , which is also Lagrangian. If one 
fixes the collection of parameters g and considers a Lagrangian manifold A+ defined 
by equations 

(x l -TT(Sl L \g)w^(P;g))\ l 

/ f) ( 1//n i 

f Vd ^-^—k B \nZ^ +d {^g) 



dx d - 



+i 



0=(yi,---,Vd) 



0. 



0. 



<x 1/n ) 



where i = 1, 2, . . . , d (i.e. the generating function is kB^Z k ^ +1 . . . , y d ); g) and 
the focal coordinates are (yi, . . . ,yd,x d+ i)), then for the manifold r~ 1 (A+) one can 
claim that 



0. 



1. 



for i = 0,1, ... ,d. These d + 1 equations do not characterize r _1 (A+) completely, 
since the dimension of the manifold is d + 1, while the sum of all these equations 
yields an identity = 0. The missing condition is the normalization: 

Now let us look at the continuous case. The analogue of the Gibbs's statistical 
operator (fl9l) is as follows: 

W ? £!>(-)^):=exp(-^ 

where a (a) is a smooth function that we use instead of {r/i}f =0 in the formula 
An analogue of the phase space R^ 4 " 1 -* is formed by pairs of functions (F(-), &(■)), 
equipped with a canonical symplectic structure u := J dz(6a)(z) A (5F)(z). An 

analogue of the manifold A+, which we denote A+ g is described as follows. For every 
L > 0, look at 



da a(a)J r ^ L \a; g) ) , 



y { e L \g) ■= {^)\TrwffM-y,g) = iy. 



Then Aj~_ is described as 



< 9 



(F(-),a(-)) F(a) 



' 5 a (a) 



, a G 



°(-)eyP(g),L>o}. 



Since we are interested in the limit e — > (i.e. the approximation S e (-) becomes the 
Dirac's delta), one can proceed as 



w£ £ K);s) 



X 



Tr 



da' exp ( — 



da a(a)S £ (a — a') ) x 



\{5 £ {a' 3 - E\ L) {g)) \+o{e)= / daexp ( - — a(o) rfto; ff ) + o(e), 

.7=1 ^ ^ B 
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where 

Tf\a-g) :=Tr{j[6 e ( aj -Sj%))) 

is the analogue of the statistical weight of a microcanonical distribution. Therefore, 
the equations describing A+ 9 in the limit e — > become more simple. Computing 
the variational derivative, for every a — (ai, a 2 , . . . , a^), we obtain: 

F(a) = L n exp ( - £<r(a))rf )(a; <?) + 0(5), 

where cr(-) varies over ^^(j?) as L varies over K>o- This motivates the following 
construction for the phenomenological Lagrangian manifold, which we denote just 
A + . We assume, that we are given a function T^ L \a), a = (ai, 02, . . . , ad), depending 
on a parameter L > 0. Set 

y (L) := {(7(0 I | rfaexp ( - ^-a(a))r( L \a) = 1}, 

for every L > 0. The Lagrangian manifold A + is described as a collection of pairs, 

^ + : = U O) I O G 3^ (L) & Va : F(a) = L"exp ( - ^a{a))^ L \a)). 

(26) 

This manifold is an analogue of r 1 (A + ). Note, that the function r^ L ^(a) should be 
linked to the phenomenological entropy S(ai, . . . , a^; V) of a nonequilibrium state 
with volume V = L n and the values of the other extensive coordinates (a 1 , a 2 , ■ ■ ■ , a^) 
via the Boltzmann's formula: 

S(a 1 ,...,a d ;V) = k B \n(cT ( < L \a)), 

where ks is the Boltzmann's constant, and c is an arbitrary constant having the same 
units of measurement as the product a\a 2 . . . (the entropy in phenomenological 
thermodynamics is defined up to an additive constant). 

At this point one can proceed as described above and consider the quasithermo- 
dynamic fluctuation theory over the points (F(-),a(-)) in terms of the asymptotic 
expansions in the inverse rescaling parameter A -1 — > 0. Without going into details, 
one can say the following. Denote the functions (F(-), er(-)) corresponding to a point 
a G A + as (F(-;a),a(-;a)). It follows from the construction of A + , that to every 
a E A + we can associate L(a) := (/ da F(a; 5)) 1//n , and a(-,a) E y( L ( a )). Assume 
that one can find a function g\(a) and a function a\(-;a) E y^ 1 nL ^\ such that 
the quantities F(a, a,X,e) defined by 

5 

'5a(a) L L J / n .) :o) 

admit asymptotic expansions as A — > 00, 



F(a,a,X,e) := ( - k Bj ^-Tr [W<£? L(5)) (a(-), g x (a)j\) 



A 1 F(a, a, A, e) ~ ^ A l Fi(a,a,e), 
1=0 

where the coefficients Fi(a,a,e), I = 0, 1, 2, ... , are smooth enough and satisfy the 
condition that, for every fixed 5 and I, Fi(—,a,e) has a weak limit as e — > 0, and 
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for the leading coefficient (corresponding to I = 0) this limit is F(-;a). Consider the 
operators of flows Ie L \a; g) corresponding to J-"j L ^(a; g): 



i=i 



where H^(g) is the second quantized Hamiltonian of the underlying multiparticle 
system, and J^ L \g) are the flow operators f[2"3j) of the observables 8 < f J \g). The 
corresponding phenomenological flows X(a; a) can be perceived in analogy with 
as follows: 



X(a; a) = lim lim < A 1 Ti 

e— s-0 A— >oo 



^>(-);^ L) (a;s) 



L=A 1 /nL(5),g=g A (5), ( 7(-)=<7A(-;5) 

Expanding the definitions of the quantities in the square brackets, we obtain: 
^i%{a{-) ]g )% L \a-g)] = Alrjexp ( - £ /daV(a')5(a' - a)) x 

d 

*fi L \g)L n J[5 E {a j -£f{g))) + o{e). 

3=1 

Performing the integration in the exponent, and taking into account, that for 5 £ A" 
we have exp(-(£(5))V(a; a)/k B ) = F(a; a) / ((L(a)) n T^ L ^\a)) , we obtain: 



Tr 



wt L B f\a(.;S);g)^\a;g) 







i=l 



(g) 



O £ . 



Therefore, the expression for the phenomenological flows reduces to 

d 



I(a; a) 



E 

i=i 



(9a,- 



F(a; a)Mj(a; 



where 



Ui(a; a) := lim lim I A x Tr<( Ji (g) 



nti^vV^Ml 

rw(a) J 



L=L x (a),g=g x (a) 



for i = 1, 2, . . . , d. Like in the beginning of this section, the evolution of the system is 
described as a curve 7 = {5*}t on A + . There exists an equation describing this curve, 
the derivation of which is totally similar to the derivation of (121]) . see [21] for details. 
Under the assumptions of the "smallness of time-memory effects" and "smallness of 
spatial gradients", one arrives at the generalized Fokker-Planck equation: 



|F(a;5.; 



i=l 



a; a }Ui{a; a ) I — 



/ <9 t .9 A F(a; a*) 



i,i=i 



0. 
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where /Cij(a; a 1 ) are some coefficients (their role is similar to the Onsager coefficients 
^i,j( at ) i n Q21j) ). and T^(a) is the function defining the limiting manifold A + (i.e. 
the phenomenological statistical weight, see ()26l)). It can be more convenient to 
write this equation in the following form: 

d d 



-F(a;5')+^ — (>(a; 5>;(a; 5 4 )) - £ I — ^(a; 5*)— lF(a; 5*) = 0, 

' i=i 4 ' i,i=i V « i/ 

(27) 



where 



Vij(a;a) :~- 



IC itj (a; a) 



r( L ( 5 ))(a) 

are the diffusion coefficients — 1, 2, . . . , d, a G A + , a G M d ), and 

„ ( (a; 5) := ^(a; 5) + £ ^| (^ r<1<5)l (< 

are the coefficients of drift (i = 1, 2, . . . , d, a G A + , a G M. d ). If we denote, for every 
i = 1,2,..., d, 

&\a) := k B ^-hxY^ L \a), 

(these are the analogues of the inverse absolute temperature and other intensive 
thermodynamic quantities in the equilibrium thermodynamics), then we obtain 

d 

Vi(a; a) = Ui(a; a) + y^P^a; a) [3f {a)) (a) , 
i=i 

for i = 1,2, ... ,d, a G A + , a = (a 1; a 2 , . . . , a d ) G M d . The equation fT271) is more 
familiar in the phenomenological physics, since in practice, £(5') is quite often fixed 
along the evolution, and the coefficients depend very little on a and 5 and can be 
determined empirically (recall, for example, the Fourier's law for the heat transfer, 
or the Fick's law for diffusion). 

The equation ( 1271) can be perceived as a "more precise" description of a physi- 
cal system, than the one given by ( 12TT) . Suppose £(5*) = const along the evolu- 
tion. Then it is more convenient to introduce f(a;a) := (L(a))^ n F(a;a), which 
is normalized as j daf(a;a) = 1. Instead of a curve {<y t }t C A + , we now have 
a curve {5*} t C A + . This function describes the density of distribution of "fluc- 
tuations" of Ei, Ei, ■ ■ ■ , E4 (recall, that A + is a Lagrangian manifold embedded 
in a phase space with extensive coordinates (xi, . . . , Xd+i) = (E\, . . . , Ed, V) and 
intensive coordinates (yi, . . . , yd+i) = (Pi, ■ ■ ■ , fld,p), the symplectic structure is 
u + = Yli=\dyi A dxi). We have a pair of consistency conditions between the two 
descriptions. The first one is that the entropy S(E U . . . ,E d ,L n ) = JEi=iA^> 
L\ry = V x l n = const, corresponding to the Lagrangian manifold A + , is linked to the 
phenomenological statistical weight r^ L ^(a), a = (ai, . . . , a d ) associated with A + , via 
the Boltzmann's formula: 

S(E 1 , ...,E d , L n ) = k B HcT^(E ll ...,E d )), 



MICROLAGRANGIAN MANIFOLDS AND FLUCTUATIONS 



21 



where fee is the Boltzmann's constant, and c > is a constant with the same units of 
measurement as (EiE 2 . . . Ed) -1 . The second one is a condition on the mathematical 
expectation: 

Ei(ct) = I daa i f(a;a t ), 



where i = 1,2, ... ,d, (note, that this equality becomes approximate, once we use 
the Onsager's equations (12"TT) and the generalized Fokker-Planck equations (j2"7]l ). If 
we wish to consider (1271) with an initial condition F{a,a t )\t=o — L n fo(a), then for 
the rest we are not restricted in a choice of /o(-) (i.e. the starting point on A + is for 
the rest arbitrary). In particular, consider the cumulants (assuming they exist): 

C N (t):= (-i& B — j / da/(a,a*)exp (-^J^MjJ 

^ i=l 

where b = (b u ...,b d ),N= (m, . . . , n d ) G Z^ , = (0, . . . , 0) G R d , N ^ 0, and 
we use a standard notation for multi-indices: (—HcbO / db) N := rij=i( — iksd / dbj) nj . 
The second consistency condition mentioned is just 



b=0 



C N (t) = lim A- 1 ^ 



A— voo 



(N,0) 



if |AT| := Yn=i n i = !> where iV = (ni,...,n rf ), and (AT, 0) := (n l7 . . . , n d , 0). On 
the other hand, this equality is not required in case where \N\ ^ 2, i.e. the cumu- 
lants C7 Ar (t) should be perceived as derivatives of the action function corresponding 
to different Lagrangian manifolds Ait "attached" to the points a 1 G A + . More 
precisely, we need only the germs of A4t of these Lagrangian manifolds (which are 
termed the microlagrangian manifolds). Note, that this construction is similar to 
the construction used in the method of canonical operator with a complex phase 

ana. 

How do M. t change with time? More generally, we would like to consider a "de- 
formation" of the linear fluctuation theory with respect to the rescaling parameter 
A — > oo, i.e. to introduce the functions C^\t) generalizing Cjv(t), A^ G Z| , N ^0. 

How does the collection {C^\t)}N change with time? The main idea which allows 
to derive these equations is as follows: let us consider the limit A -1 — > in analogy 
with the semiclassical limit h — » of quantum mechanics (where h is the small 
parameter of the transition corresponding to the Planck constant H). 

Suppose we have a classical mechanical system with n degrees of freedom q = 
(qx, q 2 , . . . , q n ) G W 1 described by a Hamiltonian H{p, q), where p = (pi,P2, ■ ■ ■ ,p n ) 
are the canonically conjugate momenta corresponding to q. Assume that if is a 
smooth function, say, from the Schwartz space <S(R?1) of rapidly decaying functions 
at infinity. A semiclassical analogue of this system is described by a Hamiltonian 
Hh{q,p) which depends on the small parameter h — >■ of semiclassical approxima- 
tion, Hh(q,p) = H(q,p) + YH°=i((~~ ih) s /s\)H ( - s \q,p), where r is a positive integer, 
and #W G S(R^t), s = 1,2 



,r . 



A semiclassical wave function ip^q) satisfies modulo 0(h°°) the Schrodinger equa- 
tion with a small parameter h in front of partial derivatives: 



ih^t(q) = H h ( y -ih-, q y h (q) + 0(h 
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where one assumes the Weyl quantization on the right-hand side, and one defines 
0(h°°) as follows: V'K?) = 0(h°°) if and only if for all k = 1, 2, 3, ... , there exists 
C\ > 0, such that sup ?gR „ |V4(<?)I ^ C l k h k . For the semiclassical wave functions 
there exists a limit 

p\q) := lim |^(g)| 2 , 

/i-»0 

(invoke the Bohr's correspondence principle). The main idea is to try to perceive 
f(a; a 1 ) = L~ n F(a, 5') corresponding to the equation ( 1271) in analogy with p t (q). 
The small parameter h — » should correspond to A -1 — > 0. Consider first the 
Wigner's quasiprobability function associated with ip^q): 

PI(Q,P) := ^ jdq' exp ( - - |)^(? + § 

where the bar denotes the complex conjugation, and (q,p) G M.^™. It follows, that 



p\q) = dp limp* h (g,p), 



where the limit under the integral on the right-hand side is taken in the weak sense. 
Note, that p l h (q,p) is real, p^((?,p) = Ph(q,p), but, in general, does not need to be 
positively defined, and satisfies only the normalization condition J dpdq p^iqjp) = 1. 
The classical limit p t {q,p), on the other hand, satisfies p t (q,p) ^ 0. Consider now iV 
semiclassical particles with an interaction between a pair of particles concentrated 
in (q,p) and {q',p') described by 

ri (—ih) s 

V h (q,p; q\p) := V(q,p; q',p') + : — V (s \q,p; q',p'), 

where r\ is a positive integer, and V and V s , s — 1, 2, 3, . . . , are elements of S{Rg H p x 
Mqp) (i.e. Swartz functions on a 2-particle phase space). To describe the kinetics of 
this system, one introduces the 1-particle density function B? Nh (q,p) (the analogue 
of Ph{q,p)) normalized on the number of particles, 

dqdpR^^q.p) = N, 

and considers a limit iV — > oo, h — > 0, in such a way that 

Nh n = x = const > 0. 

Assume that we adjust the "geometry" of the configuration space of the system (this 
can be the radius L of a n-dimensional torus, if the extracted system is confined on 
a torus), as well as the other parameters g = (gi, gi-, ■ ■ ■ , g m ) describing the "nature" 
of the system (for instance, the radius of interaction or a parameter describing the 
external field) in such a way, that for the weak limit p l {q,p) of R l N h (q,p)/N (where 
iV — > oo and Nh n = x > is fixed), one obtains an equation of the shape 

^p\q,p) + {Hlp^p; *),f?(q,p)} ~ J\A(q,P\ *) = 0, (28) 

where { — , — } denotes the canonical Poisson bracket on the 1-particle phase space, 
with coordinates q = (gi, q2, ■ ■ ■ , q n ) and momenta p = (pi,P2, ■ ■ ■ ,Pn), {Pi, qj} = 
hj ^ N) the function H [/?*] (q, p; x) is the selfconsistent Hamiltonian, and 
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J[p*](g,p; x) is the collision integral [23]. A usual way to analyse this equation is 
via an asymptotic expansion of the collision integral in x — > and by constructing 
the BBGKY chain of equations. Observe, that if we integrate ( 128|) over p, then 
this yields an expression in which the term corresponding to the Poisson bracket 
can be associated with the first sum in ( 127|) (the sum over one index), and the term 
corresponding to the collision integral can be associated with the second sum in (1271) 
(the sum over two indices). 

To make the analogy more explicit, let us look at our thermodynamic system 
and consider along with S[ (g), ■ ■ ■ , S d L \g), the flows Jx\g), • • • , Jd L \d) defined 
in (123]) (recall, that L denotes the radius of the n- dimensional torus on which we 
consider the system, and g = (#1, #2> • • • ,9m) is a collection of parameters describ- 
ing the interaction and the external field). Recall, that it is also assumed (for 
simplicity), that {S- L \g)}f =1 mutually commute, and that {<Ji (g)}f=i mutually 
commute. Look first at the analogue of (|2~T|) (we assume that the volume V = L n is 
fixed along the evolution). One needs to consider the phase space K 4d+2 with "co- 
ordinates" (xi, . . . , x d , x d+1 , x 2d , V) and "momenta" (y u . . . , y d , y d+u . . . ,y 2d ,p), 
where for s = l,2,...,d, x s corresponds to S^ L \g), x d+s corresponds to J^ d + S (g), 
y s corresponds to the "inverse temperature" conjugate to Ss L \g), and y d+s cor- 
responds to the "inverse temperature" conjugate to Js L \g)- It is convenient to 
denote x 2c i+i = V, and y 2 d+i = V (the intensive quantity "pressure over abso- 
lute temperature" associated with the volume V). The symplectic structure is 



12i=i dyi A dxi, and the evolution is described as a curve 7 = C A # 

on a Lagrangian manifold A* C M 4d+2 , where t is time. Denote Xi(a) and ?/«(«), 
i — 1, 2, . . . , 2d+l, the coordinates of a point a G A* acquired in the ambient phase 
space. The curve satisfies the following system of equations (in analogy with (12T1)): 



= Ma t) + jp LiJ (0*^(0*), (29) 

where i = 1,2,..., 2d, and Jj(a) and Ljj(a), i,j e [2d], are given functions on 
A* (the flows and the Onsager coefficients). These functions can be perceived as a 
limit A — > 00 of more complicated quantities associated with rescaled systems. In 
particular, for J s (a), s G [d], one needs to consider the flows of S\ (g) given by the 
operators ( |23i) and for J d+S (a) one needs to consider the flows of the flows Js L \g), 

&) ■= -M L \g)Jj L \g)} ={W L \g)X\n {L \g)A L \9) 



h 



h 



h 

where s G [d], (we have a family of underlying systems parametrized by L and g). 
If we denote £ { dJ ? s {g) := J s {L \g), s G [d], then the space of relevant ensembles is 
described by the Gibbs's operators 

1 2d 

g) ■= -n — ex P (-^-E ^ L) (s)) , 

where = . . . , (3 2d ), and zj^(f3; g) is determined by the normalization condition 
Triuj^ (/?; g) = 1 (assuming the trace exists). For every a G A* there exist functions 
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(3\(a) and g\(a) depending on the large rescaling parameter A — > oo in such a way, 
that 



Ji(a) = lim <^ A Tr 

A— >-oo 



i=(Ax 2d+1 (a))i/",/3=/3 A (a),p= gA (a) 



where i G [2d]. The coefficients L it j(a) can be perceived in a similar way as a limit 
of more complicated expressions (involving, for example, the projection operators of 
Mori [15] and Zwanzig [25j). It is important to point out, that the system (1291) is 
valid in the approximation of "small spatial gradients". To compute the "flows of 
the flows", one applies the commutator (i/h)\H^ L '(g), — ] twice, and, therefore, we 
should put 

J d+S (^) = 0, (30) 

where s£ [d], along the evolution curve 7 = {a t } t C A*, see [25] for details. The 
condition (I3"U|) simply restricts the possible curve 7 and implies, that the Onsager's 
matrix ||Ljj(a*)||f cannot be arbitrary. 

Now, if we wish to describe the "fluctuations" of the values of x^a*), 2/i(a;*), 
i G [2d], corresponding to a point a* G 7 C A* of the evolution curve, then we should 
introduce a "phase space" formed by the pairs (F(x),a(x)) (F is a distribution, and 
a is a function), where x = (x\, x 2 , • • • , x 2 d) G M. 2d . For the points a G A#, we have 



d ~ 

Vi{ot) = ■q^-S(x 1 , . . . ,x 2d ; V(a)) 



=Xj(a),j£[2d] 



A* : = |J [(F(-),a(-))\a(-) G ^ (L) &Vx : F{x) = L"exp ( - ^(x))r^(x)} 



where % G [2d], and S(x\, . . . , X2d\ V) is the phenomenological nonequilibrium en- 
tropy (see [18]) (note that this is a function of extensive thermodynamic coordinates 
(xi, . . . , Xd), their flows (xj+i, ■ ■ ■ , x 2 d), and the volume V). Define the phenomeno- 
logical statistical weight r^ L ^(xi, . . . ,X2d) from the Boltzmann's formula: 

S(x 1 ,...,x 2d ;L n ) = k B \n(cf (L \x 1 ,...,x 2 d)), 

where c > is a constant with the same units of measurement as (xi . . . x 2 d)~ l , and 
construct a "huge" Lagrangian manifold A# like ( I26p . 

— < 

L>0 " B 

(31) 

where x = (xi, x 2 , . . . , x 2d ), and := {a \ J dxT^ L \x)exp(—L n a(x)/kB) = 1}- 
The difference with the formula (|26|) is only that we have 2d arguments in the 
functions, in place of d. The evolution of the system is a curve 7 = {a t }t C A* 
described by a system of equations similar to ( |27l) and ( |29l) with a condition similar 
to (ED]). 

Denote (F(x;a),a(x;a)) the coordinates of a point a G A* stemming from the 
phase space. Observe, that Js L \g) = Sj^^g), s G [d], is present in the prod- 
uct of delta functions Yl^ti $( x t ~ Q (sO) corresponding to the analogue of the 
microcanonical distribution in this case. Therefore, computing Ui(x; a) for the gen- 
eralized Fokker-Planck equation (where i G [2d], x = (37, . . . , x 2 d), a G A#), wc 
obtain u s (x;a) = Xd+ S , s G [d]. The quantities Ud+ S (x;a) correspond to applying 
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the commutator (i/h)[H^ L \g), — ] twice (i.e. "flows of the flows"), and in the ap- 
proximation of "small gradients" we should put them to zero [25]. In the end, this 
yields the following equation for F(x; 5*): 



s=l i=l 1 j=l 



t& iAx ' S,) & {X ' S,) = °- 



i,j=l 

where T>^{x\ a) are the analogue of the diffusion coefficients in (|2"7j) . and f3j L \x) are 
defined as fif\x) = (k B d/dxj)\nf (L) (x), j G [2d]. The vector field on R 2 x d with the 
components Yi(x; a) := Y^jLi ^i,ji x 'i ®*) /^j 0*0, * [2d], describes a contribution 
to the total flow associated with 8\ L \g) induced by a deviation of the system from 
a thermodynamic equilibrium. Let us assume that the divergence of this vector field 
vanishes: 

2d q 

J2—Y t (x;a) = 0. (32) 
i=i OXi 

Then, for every a G A # , we can find a function <f>(x;ot) such that Yi{x\a) = 
d<&(x;a)/dxi, i G [2d], x G M. 2d . Furthermore, let us represent this derivative as 
d&(x; a)/dxi = Y^j=i Ji,jdV(x; a)/dxj, where Jij = 6ij-d~^i-d,j-, 5 is the Kronecker 
symbol. The matrix J = \\Ji,j\\fj=i satisfies J~ l = J T = — J, where (-) T denotes 
the transposed matrix. Then we have dV(x;a)/dxi = Yl"j=i(J~ 1 )i,jd^(. x ) / & x j > an d 
the function V(— ; 5) exists since Y^=i d/dxi( Y^j=i{J~ 1 )i,jd®{ x ) /d x j) = due to 
J T = — J. So one obtains: 

, _ dV{x;a) dV (x; a) 

Y s (x; a) = — , Y d+S (x; a) = , 

ox d+s dx s 

for s G [d]. Substituting this into the generalized Fokker-Planck equation for 
F(x;5*), we obtain: 

dF(x; a*) | A f dH(x; a*) dF(x; a 1 ) dH(x; jg) dF(x; a 1 ) \ _ ^ 
& ~[ V dx d+s dx s dx s dx d+s J 



where I(x,a) := Y^i j =1 ((d/dxi)V ii j(x; a)(d/dxj))F(x; a), and 

1 d 

H(x; a) := - x d+s + V(. x \ «), 



8=1 



where x = (xi, . . . , x 2d ) G M 2d , and a G A # . It is suggested to perceive H(x; a) as 
an analogue of a self- consistent Hamiltonian in mechanics (i.e. V(x; a) corresponds 
to the dressed potential), and the quantity I(x;a) as an analogue of the collision 
integral. The assumption fl32|) about the divergence of the vector field describing 
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the thermodynamic flows is basically an assumption about applicability of a self- 
consistent picture of description, which, of course, depends on a "successful" choice 
of the collection of the basis quantities {£} L \g)} d =v 

The basic idea is now as follows: let us perceive F(x; 5) in analogy with the weak 
limit of the Wigner's quasiprobability function ph{q,p) in quantum mechanics, where 
(q,p) G ^q"p, and h — > is the semiclassical parameter associated with the Planck 
constant h~. The weak limit po(q,p) = lim^Q p h (q,p) is non-negatively defined, but 
Ph{q,p) itself can be negatively defined over some region (the measure of which 
vanishes in the limit h — >■ 0, see [5])- It is suggested to perceive F(x; a) as a limit 

F(x;a) = lim F x (x;a), (34) 

A— ¥OD 

where F\(x\a) is some function (or, more generally, a distribution), depending on 
a large parameter A — > oo, and the limit is taken in the weak sense. It is quite 
remarkable, that F x (x; a) does not need to be non-negatively defined. This can 
lead to an interesting new physical effect: the thermodynamic Bell's inequalities. In 
quantum mechanics, the violation of Bell's inequalities is related to the fact that 
the probability model in quantum mechanics is different from the probability model 
in classical mechanics: the Wigner's quasiprobability can be negatively defined over 
some region of the phase space. It follows, that once we have F x (x;a) which is 
negatively defined over some region in M 2 ^, x = (x\, . . . ,X2d), we can mimic the 
Bell's inequalities and the "entangled states" in quasithermodynamics! 

Extending the analogy between mechanics and thermodynamics, one can look at 
F\(x; a) along the phenomenological evolution curve 7 = {a t } t C A*. The general 
shape of the equation for F\(x, a 1 ) should then be as follows. Equip the affine space 
M 2d , x = (X, J), formed by the points 

x = (Xi, . . . , X d , Ji, . . . , J d ), 

with a symplectic structure fl := J2 S=1 dJ s A dX s , where the first d coordinates 
(Xi, . . . , X d ) correspond to the extensive quantities (the "energies"), and the last 
d coordinates (J 1; . . . , J d ) correspond to the intensive quantities (the "inverse tem- 
peratures"). Recall, that the volume V = L n is fixed (for simplicity). Consider the 
symmetric Fock space associated with L 2 (M^): 

J=* := C © L 2 (R 2 x d ) © L 2 (M^)® S y mm 2 _ _ _ ^ 

where © symm denotes the symmetric tensor power. Let JC\(a), a G A # , be a self- 
adjoint operator on J 7 * of the form 

cfe (m) a + (x (0) )a + (^ (1) ) • • • a + (x (m) )x 

x K x m \&)a-(x {0) )a-(x {1) ) . . . a"(a; (m) ), 

where x^ -*, a;*- 1 ), . . . ^x^ vary over M M , a ± (x) are the bosonic creation and annihi- 
lation operators on J 7 * [2], 

[a _ (x), a~(x')] = = [a + (x), a + (x')], [a~(x), a + (x')] = S(x — x'), 



M 



^a(o) = 



m=0 



ml 



dx^dx^ 
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x, x' G M 2d , M is a fixed integer, x is a real parameter (introduced for convenience), 
and the operators K x m \a) are given by the commutators with respect to the Moyal 
product, 



02 (m) 



1 1 
R {m) ( (0) _ ^ J X M _ ^ J^_- 5 



where J = \\Jij\\fj =1 is the canonical symplectic matrix, H^ l \x), m — 0,1, ... , M, 
are polynomials in A -1 . 

With this setup one can mimic the constructions of quantum statistical mechanics: 
the quasithermodynamic parameter A -1 is similar to the semiclassical parameter 
h, and x is similar to the interaction parameter g. Observe nonetheless a rather 
important difference: there is a dependence of K\(a) on the point 5 of the "huge" 
Lagrangian manifold A# described as a subset of pairs (F(-),a(-)) by fj3T|) . Consider 
the second quantized analogue of the Wigner's equation for the Weyl symbol of the 
square root of density matrix, and break the symmetry in time as in 



+ ^(5*)^ = -e(R% e - r[^), (35) 

where R\ e G J 7 * is the unknown vector, and R x eq ^ G J 7 * is fixed (the value of 
lim e _s. + o R\ e corresponding to a state of thermodynamic equilibrium), and e = +0. 

Then we can perceive the solutions F(x; 5*) of the generalized Fokker-Planck equa- 
tion (1331) as a limit ( 13^1) . where 

F x (x;a t ) = \im(Ri £ ,a+(x)a-(x)Rl £ ), (36) 

and construct the analogue of BBGKY chain of equations for the higher order corre- 
lation functions. Furthermore, if there exists a generating function of fluctuations, 



2^ := J dx exp ^— uxjF\(x; a 1 ) 



where u = (ui, u 2 , ■ ■ ■ , «2d) is a parameter varying in a neighbourhood of = 
(0, 0, . . . , 0) G R M , ux := Yld=i u i x %-. then one can compute the cumulants of fluctua- 
tions Cj^(t), M G Z|o\{0}, of the extensive thermodynamic quantities . . . ,Xd) 
and their flows (xa+i, • • • , X2d) as follows: 



C${t) = ( - ik B ^-) M k B \nZ^(u; a 1 



u=0 



Note that the function Zj^ (kt; a 1 ) (if the corresponding analytic continuation ex- 
ists) restricted to Ud+i = ■ ■ ■ = u<id = yields an analogue of a partition func- 
tion in nonequilibrium statistical thermodynamics. The asymptotic expansion of 
$[^(m; 5*) := X^ks^zj^' (iu; a 1 ) in A -1 — > is an analogue of these data in nonequi- 
librium quasithermodynamics. The nonequilibrium phenomenological thermody- 
namics corresponds to the Lagrangian manifold A*. The generalized Fokker-Planck 
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equation is a phenomenological equation in this terminology, and a link between the 
"deformed" cumulants {C$(t)}M at different moments of time t is induced by ( 135]) . 

4. Thermocorpuscles 

In the semiclassical approximation of quantum theory it is a common prac- 
tice to denote the small parameter of the asymptotic expansion as h. In this 
case it is necessary to keep in mind, that the physical value of the Planck con- 
stant should not be confused with this parameter. For example, one may denote 
fophys — 6.6262 x 10~ 27 erg s and after that write h — > 0. In quasithermodynam- 
ics, the small parameter is A -1 , where A is the rescaling parameter. In can be 
convenient to redenote it as k B since it stands in the same place as the Boltz- 
mann's constant in the exponent linking the free energy and the partition func- 
tion. If we redenote the physical value of the Boltzmann constant, for example, as 
(kB) P hys = 1-3807 x 10~ 16 erg K -1 , then we may write k B = A -1 — >■ 0. 

Intuitively, when the Planck's and Boltzmann's constants are introduced into 
physics, h p hys is related to the quantization of "properties" (for example, the energy 
spectrum), and {k B ) phys is related to the quantization of "substance". In this sense, 
the formula for the quantization of energy E of a 1-dimensional harmonic oscillator 
is similar to the "quantization" of the number of moles u, 

E — E = u h ph y S n, v = R~ x {k B ) phys N, 

where n G Z^ (the number of quanta), N G Zj> (the number of particles), E is 
the ground level of energy, u > and R > are parameters (the frequency of the 
oscillator and the universal gas constant, respectively). 

In the previous section we have introduced the creation and annihilation opera- 
tors a ± (x), x = (X, J), where X = (X 1; . . . ,X d ) are the extensive thermodynamic 
quantities (like the values of internal energies or the numbers of moles in different 
parts the system), and J = ( J 1; . . . , J^) are the associated flows. The volume of the 
system V is fixed (otherwise, we need to add one more coordinate X2d+i = V). It is 
natural to perceive these operators as the operators of creation and annihilation op- 
erators of "thermoparticles" (or, another name, could be "thermocorpuscles") . We 
work over the space Rjj rather than K.^-, or M.j, in order to avoid the discussion 
about the statistics of these thermodynamic particles. Are they thermo-fermions, 
thermo-bosons, or, perhaps, something else? Leaving this for another paper, let us 
look at two other effects. 

Effect no.l: The thermodynamic Bell's inequalities. In ( 133]) we have F(x; 5*), x = 
(X, J), which is an analogue (up to a normalization factor L n ) of the classical prob- 
ability distribution (recall, that 7 = {5*} C A* defines a curve of evolution of the 
system). In particular, for f(x; a) := L~ n F(x; a), we have J dx f(x; a 1 ) = 1, and for 
any test function <f(x) G C^°(M^ d ), such that (p(x) > 0, we have J dx cp(x)f(x; a 1 ) ^ 

0. On the other hand, for the "deformed" distribution f\(x;a) := L~ n F\(x;a), 
one does not impose the latter condition, and there can exist a test function ip\ G 
C^°(I^ d ), such that <p\(x) > 0, but f dx ip\(x)F x (x; a f ) < 0. In mechanics, the fact 
that the Wigner's quasiprobability function can be negatively defined over some re- 
gion of the phase space is responsible for the Bell's inequalities [Ij. In particular, one 



MICROLAGRANGIAN MANIFOLDS AND FLUCTUATIONS 



29 



may introduce a thermodynamic wave function, considering, for example, a Cauchy 
problem with an initial condition stemming from an element ip\(X) G L 2 (M^-), 

/(*;S')U = p^rip JjX' exp ( - ^JX')i;(x - + f ), 

where X = (xi, . . . , x^), J = ( Ji, . . . , Jd)- The parameter A -1 in this formula is an 
analogue of the semiclassical parameter h in quantum mechanics. If we take, for 
instance, a WKB-type function 

^ A (X)=exp(^S(X))^(X), 

where S(X) is real, and <p\(X) admits an asymptotic expansion in the powers 
of the small parameter A -1 , then we should interpret the derivatives dS(X)/dXi, 
i = 1,2, ... ,d, as the flows J = (Ji, J 2 , ■ ■ ■ , Jd) induced by X = X 2 , . . . , X d ). 
The limit A -1 — > corresponds to the nonequilibrium thermodynamic Hamilton- 
Jacobi equation. In other words, the phenomenological thermodynamic flows are 
the analogues of classical mechanical momenta. Intuitively, the leading coefficient 
liniA^oo fx{X) is something similar to a Gaussian exponent concentrated near a 
point. Observe that the quadratic function corresponding to the power of this ex- 
ponent does not contain the small parameter A -1 . Considering a thermodynamic 
system built from several similar subsystems, one may construct an entangled state 
and violate the Bell's inequalities. The "observables" in this case are represented 
by self-adjoint operators on L 2 (M^) depending on a parameter a G A*, and the 
inequalities mentioned yield a condition on the fluctuations of the measured values 
of Xi and Jj, i,j = 1,2, ... ,d. 

Effect no. 2: The "deformed" Boltzmann's H -theorem. The effect no.l discussed 
above corresponds to the first quantization of phenomenological nonequilibrium ther- 
modynamics, where the Boltzmann's constant ks plays, in a certain sense, a role 
similar to the Planck's constant ft. In statistical mechanics the particles interact with 
each other (for example, via an interaction potential). The parameter g in front of 
the interaction potential is sometimes termed the external Planck constant g = h~ ext , 
since the commutation relation [y/ga + (q), y/ga~{q')\ = g5(q — q') for the bosonic 
creation-annihilation operators in configuration space points q, q' G W 1 are similar 
to [bf, bj] = HSij, where bf = (qj =F fd / dqj) / \/2, i,j = 1,2, ... ,n. In the previous 
section we have seen, that it is natural to introduce an interaction between the ther- 
mocorpuscles (the collision integral). Let F\(x; 5*), x — (X, J), be of the shape fl36l) 
corresponding in the limit A — > oo to a solution of the generalized Fokker-Planck 
equation f )33|) . One perceives F x (x; a 1 ) as an analogue of a one-particle kinetic func- 
tion in quantum statistical mechanics. If F\{x;a t ) corresponds to a symbol of a 
A _1 -pseudodifferential operator, 

F^^F^X^iX- 1 — ;^, 
where one uses the Weyl quantization, then it is possible to consider a quantity 
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In analogy with the Boltzmanns theorem, the presence of the collision integral im- 
plies, that 

along the nonequilibrium evolution curve 7 = {a t }t C A*. This inequality should 
hold not just in the limit A — > 00, but for all values of A. 

As a final remark, it is of interest to point out the following. A collection of 
extensive thermodynamic quantities X = (X\,X2, . . . ,Xd) is normally a collection 
of local internal energies or local densities of chemical substances. If we consider 
the system together with its environment (i.e. the thermostat) as one big system, 
then we have conservation laws (for example, the total energy is fixed). It follows, 
that instead of the canonical symplectic structure u = Yli=i dJi A dX{ on the whole 
phase space ifj, one may reduce the problem to a smaller phase space separating 
the integrals of motion (in mechanics this step is termed a reduction of dynamics to 
a coisotropic sub manifold) . The reduced symplectic structure can be perceived as a 
kind of "thermodynamic gauge field" created by the environment. 

5. Discussion 

The main motivation of the present paper is a striking (though underappreci- 
ated) analogy between the roles played by the Planck's and the Boltzmann's con- 
stants in the mathematical formalism of theoretical physics. These two constants 
were introduced into science at more or less the same time (the beginning of the 
XX-th century) and in a sense capture the Zeitgeist (the spirit of the age). The 
concept of a thermocorpuscle, i.e. a thermodynamic particle with the thermody- 
namic forces X = (X 1? X 2 , . . . , Xj) in place of coordinates (internal energy, number 
of moles, volume, etc.), and the associated nonequilibrium thermodynamic flows 
J = (Ji, J 2 , . . . , Jd) in place of momenta, is, in a sense, quite natural if one investi- 
gates the analogy between the Heisenberg's uncertainty principle and the Einstein's 
quasit her mo dynamic fluctuation theory. Informally speaking, the Planck's constant 
h and the Boltzmann's constant ks are very similar concepts, but corresponding to 
different "hierarchical" levels of organization of matter. 

The symmetry between the mechanical and thermodynamic pictures is very much 
in line with the general philosophy of I. Prigogine [18] who tried to introduce the 
"irreversibility on a microscopic level of description" . A thermocorpuscle is a coun- 
terpart of a quantum particle, with h replaced by k&. It is of interest to point out, 
that the analogy between h and ks is a general motivation behind the mathemat- 
ical constructions in [21] and [22], where one describes a kind of "noncommutative 
neighbourhood" around the semiclassical and quasithermodynamic parameters. Per- 
haps even more complicated new theoretical structures are needed to capture the 
nonequilibrium statistical physics in all its aspects. 

At this point it is possible to claim that the "quantization" of phenomenological 
thermodynamics (with &b i n place of h) and the "second quantization" of the quan- 
tized thermodynamics appear to be quite reasonable constructions. In particular, 
conceptually it is important to distinguish between the phenomenological thermo- 
dynamics, statistical thermodynamics, and the thermodynamics that lives between 
these two theories, the quasithermodynamics. This is similar to mechanics: there 
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is classical mechanics, quantum mechanics, and the mechanics existing between the 
two theories, the semiclassics [11]. If we look at the deformation quantization the- 
ory [7], then it is a standard practice to denote the "formal" parameter in the star 
product as h. Now one can see, that ks is just as good: 

oo 

f(x, J) * kB g(x, J) = f(x, j) g (x, J) + J2 k B B sU, g)(x, J), 

8=1 

where f,gE (7°°(R^ j), B s are bilinear differential operators, s = 1, 2, 3, . . . . 

In our days, due to the progress in the field of nanotechnologies, one can directly 
work with systems consisting of relatively small numbers of quantum particles (con- 
sider microtransistors in electronics, molecular machines in biological physics, etc.). 
For example, the size of a single transistor on a microchip is measured in hundreds 
of atoms, but certainly not in the scale of the Avogadro number. The usual empiri- 
cal laws stemming from phenomenological thermodynamics, for example, the linear 
Ohm's law for the electric current, break down at this scale. On the other hand, the 
systems are still too big to be peceived purely quantum mechanically. One needs a 
kind of "nanothermodynamics" in this case. If we look at mechanics and imagine 
that we go from a classical description to a quantum description via the semiclassics, 
then one can say that we first construct a Schrodinger equation with a semiclassical 
parameter h — > in front of partial derivatives, and then formally substitute h = 1 
(this corresponds to the "real" Schrodinger equation with h). It looks like a natural 
speculation, that if one performs a similar step in quasithermodynamics (substitute 
A" 1 = 1), then one obtains the "real" quantized thermodynamics. 
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